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Self-consistent system of nonlinear spinor field and Bianchi I (BI) gravitational one with time dependent gravitational 
constant (G) and cosmological constant (A) has been studied. The initial and the asymptotic behaviors of the field 
functions and the metric one have been thoroughly investigated. Given A = A /r 2 , with r = yf—g, G has been 
estimated as a function of r. The role of perfect fluid at the initial state of expansion and asymptotical isotropization 
. . . process of the initailly anisotropic universe has been elucidated. 

o 

^ ' 1. INTRODUCTION 

7T" 1 ! Einstein's theory of gravity contains two parameters, considered as fundamental constants: Newton's gravitational 

ry\ • constant G and the cosmological constant A jjj . A possible time variation of G has been suggested by Dirac M and 

extensively discussed in literature &-R]- The "cosmological constant" A as a function of time was studied by many 

authors. Chen and Wu M advocated the possibility that the cosmological constant varies in time as 1/i? 2 , with R 

being the scale factor of Robertson- Walker model. Further Abdel- Rahman considered a model with the same kind 

^. , of variation, while Berman et al. MT12J stressed that the relation R ex t plays an important role in cosmology. 

£S) ' Berman and Gomide [|13| also showed that all the phases of the universe, i.e., radiation, inflation, and pressure- free, 

(^> may be considered as particular cases of the decelaration parameter q = constant type, where 

o. 

ON ■ q = -RR/R 2 , (1.1) 

o ; 

^D . where dots stand for time derivative. This definition was extended by Singh and Agrawal E3 to the Bianchi cosmo- 
logical models. Recently we studied the behavior of self-consistent nonlinear spinor field (NLSF) in a Bianchi I (B-I) 
Q \ universe Eq] that was followed by the study of the self-consistent system of interacting spinor and scalar fields Em . 
Hi These studies were further extended to more general NLSF in presence of perfect fluid 0,|8| ■ 

; I | The aim of this paper is to extend our study with time dependent gravitational constant G and cosmological 
W)' constant A in Einstein's equation. 



- * — i 



2. FUNDAMENTAL EQUATIONS AND GENERAL SOLUTIONS 

The Dirac spinor field is given by the Lagrangian 



L= l - 
2 



^v^v - V M V7^ 



impip + L N . (2.1) 



The nonlinear term Lm describes the self-interaction of a spinor field and can be presented as some arbitrary functions 
of invariants generated from the real bilinear forms of a spinor field having the form 

S = tpip, P = #7 5 V>, w M = (VV'i/O, Atl = (V'7VV'), T^ = (i/ja^ip), 
where a^ v — (i/2)[7 M 7 l/ — r f^\. Invariants, corresponding to the bilnear forms, look 

I A = A^A» = (^ 7 VVO fl^$7VV0> !t = V T"" = (^^) g m g vP (^G^i,). 

According to the Pauli-Fierz theorem p9| , among the five invariants only / and J are independent as all other can 
be expressed by them: I v = —Ia = I + J and It = I — J- Therefore we choose the nonlinear term Lm = XF(I, J), 
thus claiming that it describes the nonlinearity in the most general of its form. Here A is the coupling constant. 
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The NLS F eq uations and components of the energy-momentum tensor for the spinor field corresponding to the 
Lagrangian (2.1) are 



and 



where 



i-f^Vfiip - $V + iGj 5 ip = 0, 
iV^^ + §4> - it/^7 5 = 0, 



T t = -,9 PV in^^ + h^^ ~ V^ 7 vV» - V v fy^ - 8?L sp + T p a 



~ fi (ra) ' 



(2.2a) 
(2.2b) 



(2.3) 



$ 



V 



-2XS 



dF 



g = 2\p 



dF 

dJ' 



Here T p , , is the energy-momentum tensor of a perfect fluid. For a Universe filled with perfect fluid, in the 



~/i(m) 



concomitant system of reference (u° = 1, u l — 0, i — 1, 2, 3) we have 

T Mm) = (P + e)%«" - <%p = (e, -p, -p, -p), 



(2.4) 



where energy e is related to the pressure p by the equation of state p = £ e. The general solution has been derived 
by Jacobs |20[ . Here £ varies between the interval < £ < 1, whereas £ = describes the dust Universe, £ = | 
prese nts radiation Universe, | < £ < 1 ascribes hard Universe and £ = 1 corresponds to the stiff matter. In (2.2) 
and (2__3) V M denotes the covariant derivative of spinor, having the form |2l[ | 



V ^ = S" r ^' (2 ' 5) 

where r /i (a;) are spinor affinc connection matrices. 

Einstein's field equations with variable cosmological and gravitational "constants" A and G are given by 

K-\ SSR = -^G{t)T£ + A(t)% (2.6) 

where B% is t he Ricci tensor; i? = g^ v Rg^is the Ricci scalar; and T^ is the energy-momentum tensor of matter field 
given by (|2.3|). From the divergence of (|2.6j) we get 



8ttG ^ + 87rG(T^) - A, p^ = 0, 



(2.7) 



The Bianchi I model is given by 



ds 2 = dt 2 - a 2 (t) dx 2 - b 2 (t) dy 2 - c 2 {t) dz 2 . (2.8) 

We stu dy t he space-independent spinor fields, h enc e T$ is the function of t alone. Taking this into account for the 
metric ( [2.8| ), the Einstein's equations (Z.f) and (2/7) reduces to 



be 



IttGT" 



SttG 



b c be 
c a ca 

c a ca 

a b a b 

a b a b 

a b be ca 

-T + T-+-- 

a b be ca 

!? + <- + £ + 
\a b 



ZirGTl - A, 



- + -- =87tGT 2 2 -A, 



= 87rGT 3 3 


-A, 






= 8ttGT 


-A, 






c) a 


+ T 2 


h 

b 


, T 3 C " 



-A = 0, 



(2.9a) 
(2.9b) 

(2.9c) 
(2.9d) 

(2.10) 



where poin ts denote differentiation with respect to t. If we suppose the energy conservation law T^ = to hold, 
then ( 2.10J ) reduces to 
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b c\ „i a m9 b 






= 0, 



(2.11a) 
(2.11b) 



8ttGT ° - A = 0, 

Let us now go back to the spinor field equations ( |2.2| ). Using the equalities P,|22|] 

9ny{x) = e*(x)et(x)r) ab , ^(x) = e°(a;)7 Q , 

where r^ab = diag(l, — 1, — 1, — 1), j a are the Dirac matrices of Minkowski space and e"(x) are the set of tetradic 
four- vectors, we obtain the Dirac matrices 7^ (x) of B-I space-time 

7° =7°, 7 1 =7 1 /«W, 7 2 -7 2 /&W, 7 3 =7 3 Mt), 



70 = 70, 7i = iMt), 72 = 72&(*): 73 = 73c(i). 



The r^(x) matrices are defined by the equality 



which gives 



T,(x) = -g pa (x){d,e b s e p b -T^)r 1 s , 



r = o, r x = la(t)tf, T 2 = \b{t)ff, T 3 = h(t)ff. 



(2.12) 



Flat space-time matrices we choose in the form Caj 

10 
10 
0-10 
0-1 

/ 0-t 

_ 2 [ i 

7 ~ i 

\ -i 

Defining 7 s as follows, 



1 

10 

0-100 

-10 

10 
0-1 
-10 
10 



we obtain 



Defining 



-^E^ V(jp Yl v l' 7 l p , E^p = \T r ge pv(jp , £0123 = 1, 



7 = — i\/ — .97 111 



-if-1-2-3 = - 5j 



/ 0-10 
0-1 
-10 

V -1 



r(t) = a(t)b(t)c(t) = y/=g, 
we rewrite the equation ( 2.2a) together with (2.5) and ([2.12) 



^°(~di + 27 )V'-^ + ^7 5 V'= ( ' 



(2.13) 



(2.14) 



Defining E/j(u) = y/ripj(t), where j = 1,2,3,4, and v — J Qdt from (2.14) one deduces the following system of 
equations 
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ui + i($/g)u 1 -u 3 = o, 

[7 3 - i($/G)U a + Ui = o, 

U' 2 + i($/G)U 2 - U A = 0, 
Ui - i{$/Q)Ui + U 2 = 0, 

where prime denotes differentiation with respect to v. 

From ( p.2[ ) one can write the equations for 5 = ^ptfj, P — i-ip^ip and A — iprf^ip 

S Q -2gA Q =0, 

P -2$A = 0, 

A + 2<PP + 2gS Q = 0, 

where Sq = r5, Pq = 7"P, and Aq = tA, leading to the following relation 

S 2 + P 2 + A 2 = C 2 /t 2 , C 2 = const. 



(2.15a) 
(2.15b) 



(2.16a) 
(2.16b) 
(2.16c) 

(2.17) 



Let us now solve the Einstein's equations. To do it we first write the expressions for the components of the energy- 
momentum tensor explicitly. Using the property of flat space-time Dirac matrices and the explicit form of covariant 
derivative V M , one can easily find 



T ° =mS- XF(I, J) + e, Tf = if = T 3 3 = VS + QP - XF(I, J) - p. 
Summation of Einstein equations ( [2. 9a ), (2. 9b), (2. 9c) and (2.9d) multiplied by 3 gives 

- = 12ttG(t^ + T °) - 3A = 12ttG (mS + VS + QP- 2XF(I, J) + e - p) - 3A. 



(2.18) 



(2.19) 



Let us express a,b,c through r. For this we notice that subtraction of Einstein equations ( [2.9b[ ) and ( 2.9a ) leads 
to the equation 



a b he be d I a b 
a b ac be dt\a b 



a b\ fa b c 
a b J \ a b c 



0. 



with the solution 



Analogically, one finds 



a ( ' f dt\ 

— = Diexpl X\ I — I , D\ = const., X\ — const. 



^ = D 2 cJx 2 f^), b - c =D^JxJ d ± 



(2.20) 



(2.21) 



(2.22) 



where D 2 , D 3 , X 2 , X 3 are integration constants. In view of (2.13) we find the following functional dependence between 
the constants D\, D 2 , D 3 , X\, X 2 , X 3 : 

D 2 = D 1 D 3 , X 2 =X 1 +X 3 . 

Finally, from (2.21) and ( 2.22J ) we write a(t),b(t), and c(t) in the explicit form 

"2Xi+AT 3 f dt 



a{t) = ( J D2 jD3 )i/3 T i/3 exp 
b{t) = ( J D- 1 J D 3 ) 1 / 3 r 1 / 3 exp 



c(t) = (DiDly^T^exp 



3 J r(t) 

X!-X 3 f _dt_ 

3 J r(t) 
Xi + 2A 3 f dt 



(2.23a) 
(2.23b) 
(2.23c) 



Thus t he system of Eins tein's equat ions is complet ely in tegrated. Let us now go back to the equations ( 2.11a| ) and 
( [2.11b ). In view of ( 2.18| ) and (2.13) we find from ( [2. 11a ) we find 
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e+(e+p)- + ^So-gP o = 0. 

T 



5 
(2.24) 



On the other hand from (2.16a) and (2.16b) we have 

*5 - gPo = 0. 



Taking this into account and also the equation of state p = C,e, < C < 1 from (2.24) we find 



e = 



£q 

rl + C^ 



Ceo 

T !+C : 



(2.25) 



where £o is the integration constant. Let us now define G. Taking into account that G = rdG/dr and A = fdA/dr 
we rewrite (2.11b) as 



~ T odG dA 



On the other hand, inserting a, b 7 c from ( 2.25 ) into ( 2.9d ) we obtain 

r 2 X 
87rT n °G = —z - —s + A. 



3r 2 3r 2 



where X = X{ + X-lX 3 + X§. Dividing ( |2.26| ) by ( |2.27D we find the following equation for G 

dG/dr _ 3t 2 <9A/<9t 



G 



f 2 -A- + 3r 2 A' 



(2.26) 



(2.27) 



(2.28) 



Now, A is a given function of r, namely, A = Ao/r 2 as well as T^ and T§. Then (2.19), multiplied by 2f can be 
written as 



2ff = [2(l2irG(Tl +T§) -3A)r]f = ^(t)t = ^{t) 



Solution to the equation ( 2.2E ) we write in quadrature 

cIt 



vni 



(2.29) 



(2.30) 



Giving the explicit form of F(I, J), from (2.30) one finds concrete function r(t). Once the value of t is obtained, one 
can get expressions for components ipj(t), j — 1, 2, 3, 4. 

In what follows, we analyze the solutions obtained previously. In J17| we gave a detailed analysis of the problem 
for different F(I, J). Here we give a brief a ccoun t of that. 

Setting F = F(I), i.e. when Q = from ( 2.16a ) one finds 



S = — , C = const. 

r 



For the spinor field in this case we obtain 



A(t) = (C r /V^) 



-%Cl 



r = l,2, Mt) = (Ci/V^)e m , 1 = 3,4 



where fi = j $di and C r and C; are integration constants such that Cq = C\ + C\ — C\ — C\ 
If one sets F — F(J) and puts m = 0, i.e. when $ = from (2.16b) one finds 



P(t) 



Do 



D Q = const. 



In this case for spinor field we obtain 



V-i = (l/V^)(Die lv + iD 3 e~ m ), ^ 2 = {l/^){D 2 e m + iD 4 e~ iv ), 
V>3 = {l/y/ : F){iD 1 e iv + D 3 e~ iv ), i> 4 = {\/^){iD 2 e m + D 4 e~ m ). 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



where D Q = 2 (D\ + D\ - D\ - D\). 
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Let us note that, in the unified nonlinear spinor theory of Heisenberg, the massive term remains absent, and 
according to Heisenberg, the particle mass should be obtained as a result of quantization of spinor prematter Q . 
In the nonlinear generalization of classical field equations, the massive term does not possess the significance that it 
possesses in the linear one, as it by no means defines total energy (or mass) of the nonlinear field system. So our 
consideration massless spinor field is justified. 

Another choice of nonlinear term is F — F(K±), K+ = I + J = I v = —Ia, K- = I — J = It- In the case of 
massless NLSF one finds 

S 2 ± p2 = D±^ (235) 

T 

In all the cases mentioned above we mainly found t — at for small t guarantying anisotropic behavior of the 
universe at initial state, while r = fit 2 as t — * oo which is in accord with present day isotropic state. Note that this 
result was obtained for G constant. Here a and (3 are constants. As one sees, for r = at as t — > 0, the solutions of 
spinor field are initially singular. But for some special cases, it is possible to obtain the solutions those are initially 
regular pj,[l8|, but it violates the dominant energy condition in the Hawking-Penrose theorem p5| . Note that one 
comes to the analogical conclusion choosing F = F(K), K — I. J. 



Now, setting A = A /t 2 from ( |2.28| ) we find 

G = C/T 6Ao/(a2 ~ x+3Ao \ r = at, C = const., 

and 

/ ABt A 6Ao/(*-3A ) 

G = D[ - — , r = Bt 2 , D = const.. 

If w e cons ider A = Ao/r 2 and G— constant, then the conservation law T£ — does n't ho ld separately, as in that 
case ( 2.11b| ) leads to A = const., which contradicts our assumption. In this case from ( |2.10 ) we find 



with the solution 



i+(l + T -=-^f (2.36) 



* 2A ° l row 

e=— V2 (2.37) 



Setting F — K n with K = {I, J, (/ ± J), IJ} from (2.3C) we conclude that even in presence of time dependent A 



in the Einstein's equation perfect fluid plays no role at the early stage of expansion as well as isotropization of BI 
universe leaving it to the nonlinear spinor term in ( |2.l| ) which confirms our claim made in p^7|jlq | . 

Finally, we see what heppens with the system in absence of spinor field. As one sees, in this case the relation 



£ = Eq/t 1 ^ takes place. Given A = Aq/t 2 from (2.11b) for G one finds 



G 4ne (l- Qr^O (2 ' 38) 



3. CONCLUSIONS 



Exact solutions to the NLSF equations have been obtained for the nonlinear terms being arbitrary functions of 
the invariant I = S 2 and J — P 2 , where S = ■ijji/j and P — iip^f 5 ip are the real bilinear forms of spinor field, for B-I 
space-time. It has been shown that introduction of time dependent A term in Einstein's equation and consideration 
of gravitational constant to be a function of time do not effect the initial singularity and asymptotic isotropization 
process which is dominated by the nonlinear spinor term in the Lagrangian. It has also been shown that the results 
remain unchanged even in the case when the B-I space-time is filled with perfect fluid. 
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